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Abstract Optical techniques are finding widespread 
use in analytical chemistry for chemical and bio-chemical 
analysis. During the past decade, there has been an in- 
creasing emphasis on miniaturization of chemical anal- 
ysis systems and naturally this has stimulated a large 
effort in integrating microfluidics and optics in lab-on- 
a-chip microsystems. This development is partly defin- 
ing the emerging field of optofluidies. Scaling analysis 
and experiments have demonstrated the advantage of 
micro-scale devices over their macroscopic counterparts 
for a number of chemical applications. However, from 
an optical point of view, miniaturized devices suffer dra- 
matically from the reduced optical path compared to 
macroscale experiments, e.g. in a cuvette. Obviously, 
the reduced optical path complicates the application of 
optical techniques in lab-on-a-chip systems. In this pa- 
per we theoretically discuss how a strongly dispersive 
photonic crystal environment may be used to enhance 
the light-matter interactions, thus potentially compen- 
sating for the reduced optical path in lab-on-a-chip sys- 
tems. Combining electromagnetic perturbation theory 
with full-wave electromagnetic simulations we address 
the prospects for achieving slow-light enhancement of 
Beer-Lambert-Bouguer absorption, photonic band-gap 
based refraetometry, and high-Q cavity sensing. 



1 Introduction 

Optofluidies is a new branch within photonics which 
attempts to unify concepts from optics and microflu- 
idics [212]. The purpose of having these two fields play 
in concert is at least twofold; 

i) Optofluidies may provide new means for doing optics, 
including fluidic tuning of optics 



ii) Optofluidies may find widespread applications in lab- 
on-a-chip systems [71I51I51ITU] and analytical chemistry 
in general, where optical techniques are being fre- 
quently used for chemical and bio-chemical analy- 
sis [H]. 

The increasing awareness from the lab-on-a-chip com- 
munity of optofluidies is to a large extend driven by the 
past decade's emphasis on miniaturization of chemical 
analysis systems |12| which has naturally stimulated a 
large effort in integrating microfluidics [T^ITl] and op- 
ties in lab-on-a-chip microsystems [71[5]. While a broad 
variety of phenomena, properties, and applications ben- 
efit from the down-sealing in size |12j , the conditions for 
optical sensing and detection are impeded tremendously 
by the reduced optical path length in lab-on-a-chip sys- 
tems. Panel (b) in Fig. [T] illustrates a typical lab-on-a- 
chip implementation of an absorbanee cell. The optical 
path length L is often reduced by several orders of mag- 
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Fig. 1 Panel (a) illustrates a typical macroscopic cuvette 
while panel (b) shows a microscope image (top-view) of 
an equivalent lab-on-a-chip implementation of a microflu- 
idic channel (vertical direction) integrated with planar op- 
tical waveguides (horizontal direction). Courtesy of K. B. 
Mogensen and J. P. Kutter (MIC - Department of Mi- 
cro and Nanotechnology, Technical University of Denmark, 
www.mic.dtu.dk / microTAS) . 
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nitudc compared to typical macroscopic counterparts, 
such as the cuvette shown in panel (a). A typical size re- 
duction by two orders of magnitude will penalize the op- 
tical sensitivity in an inversely proportional manner, e.g. 
in a quantitative concentration measurement the lowest 
detectable chemical concentration will be increased by 
the same order of magnitude [15| . 

The above problem may conveniently be explained by 
the reduced light-matter interaction time r given by the 
Wigner-Smith delay time. For a homogeneous system of 
length L we have 

T ^ L/C, (1) 

with c being the speed of light. For small systems the 
photons will thus eventually have too little time for in- 
teracting with the chemical species and only very high 
concentrations can be quantified. Many applications in- 
volve i) minute sample volumes containing ii) very low 
concentrations of e.g. bio-molecules. There is thus obvi- 
ously a strong call for miniaturized systems acknowledg- 
ing the former without jeopardizing optical approaches 
quantifying the latter. In this work we explore the use 
of photonic crystal concepts for such a purpose. 

Photonic crystals [TBlfTT] constitute a class of ar- 
tificial electromagnetic structures offering highly engi- 
neered dispersion properties not offered to us by any 
elements of Nature herself. Photonic crystals are porous 
structures, but with the voids arranged in a highly reg- 
ular way, see Fig. [5]). This makes electromagnetic radia- 
tion interact strongly with the matter for wavelengths A 
comparable to the periodicity A of the photonic crystal. 
Photonic crystals are strongly dispersive environments 
supporting a number of phenomena including photonic 
band gaps and slow- light propagation [T81IT9]. As we 
shall see, liquid-infiltrated photonic crystals will over- 
all inherit the unusual dispersion properties from their 
non-infiltrated counterparts, thus also changing the way 
light interacts with bio-molecules dissolved in the liq- 
uid. The effects can be quite spectacular compared to 
light-matter interactions in a spatially homogeneous liq- 
uid environment. 

Previous experimental work on liquid-infiltrated pho- 
tonic crystals may roughly be classified as studies em- 
phasizing either fluidic tuning of optical properties [3l 
I20j or work focusing on optical sensing in fluidic en- 
vironments [TUpril^l^l^ipSpSPTll^ . On the theo- 
retical side there has also recently been a considerable 
attention to tuning [29] and sensing aspects [30 l l3T | [32 t 
[33 l [34 l l35 l [36] In passing, we note that somewhat paral- 
lel to the above development there has been an effort in 
fluidic tuning of photonic crystal fibers [57[I551[5^ and 
microstructured fibers have recently been integrated in 
lab-on-a-chip systems for sensing applications [40j . Fi- 
nally, infiltrating with liquid crystals [41' seems to be an 
interesting direction for dynamically tunable photonic 
crystals |42 and photonic crystal fiber devices [33]. 

In this paper we outline our recent theoretical and 
numerical work [33ll341l35[l36| on liquid-infiltrated pho- 




Fig. 2 Schematic of the dielectric function variations in a 
photonic crystal. High-index material with dielectric function 
£d is shown in gray while the liquid with dielectric function 
£i is indicated by blue coloring. The dashed curve indicates 
the dielectric function variation in the absence of liquid. 

tonic crystals including potential applications for en- 
hanced light-matter interactions in lab-on-a-chip systems. 
The remaining parts of the paper are organized as fol- 
lows: in Sec. [2] we introduce our theoretical formalism, in 
Sec. [3] we address the dielectric properties of physiologi- 
cal liquids, in Seed] we discuss general aspects of modes 
in liquid-infiltrated photonic crystals, and in Sec. [5] we 
give examples of potential applications in optical sens- 
ing. Finally, in Sec. [6] a discussion and conclusions are 
given. 

2 Formalism 

We consider photonic crystals made from a single high- 
index material (with dielectric function Sd) which can be 
of either void-like or pillar-like nature. By infiltrating the 
photonic crystal by a liquid (with dielectric function ei ~ 
nf) we have a spatially dependent dielectric function 

e{r) = I (2) 

where and "fi are the domains of dielectric material 
and liquid, respectively Fig. [2] illustrates an example of 
a planar photonic crystal and also indicates how the di- 
electric function changes upon liquid infiltration. 

2.1 The wave equation 

For the dielectric problem constituted by Eq. ([2]), the 
Maxwell equations lead to the following generalized eigen- 
value problem for an electrical field with a harmonic time 
dependence [T8llT9| 

VxVx ]£;,„)= e^|£;„) (3) 

with cj„i being the angular eigenfrequency of the mth 
eigenfield \Em'). We have for simplicity suppressed the 
spatial dependence of the dielectric function e(r). Since 
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the eigenvalue problem is of the generalized form, the 
eigenficlds are orthonormal in the sense that 



(-E„|e|-E„i) — S„ 



(4) 



where we have used the Dirac bra-ket notation with inner 
products between functions g and h defined according to 



dr g*{r)h{r) 



(5) 



with the integral extending over all space unless other- 
wise indicated. 



2. 2 Block modes of photonic crystals 

For a photonic crystal the dielectric function e of Eq. ([2]) 
is periodic in space, i.e. e(r + fi) = e(r) with R being a 
lattice vector, see Fig. [2] The solutions to Eq. ^ arc in 
this case known as Bloch modes [T8l[T9] and are charac- 
terized by both a mode index m and a Bloch wave vector 
K. Introducing the explicit time and space dependencies 
we have 

|-E'„i,K(r, t)) = exp [z (k • r - ujmt)] \E„,^^{r)) (6) 

where we have introduced the Bloch function | -Em 
which inherits the periodicity of e(r). Typically, one has 
to rely on numerical approaches to the calculation of the 
Bloch modes of Eq. ([3]) in a photonic crystal. For this 
purpose we in this work employ a plane-wave expansion 
method with periodic boundary conditions pi] yielding 
the dispersion relation u}m{i<-)- In general, the dispersion 
in a photonic crystal differs significantly from the free- 
space dispersion relation uj{k) = ck. 



2.3 Perturbation theory 

Many aspects of sensing will be adequately addressed 
with the aid of perturbation theory and in the following 
we briefly outline the basic equations of first-order elec- 
tromagnetic perturbation theory. We assume we make 
a small frequency-independent perturbation Ae to the 
dielectric function, i.e. 



Ae 



(7a) 



The perturbation is small in the sense that Ae <^ e 
and we would like to know how this is going to change 
the cigenfrcqucncies and eigenficlds. Expanding the fre- 
quency and the field in a similar manner 



\e) \e) + \AE) 



(7b) 
(7c) 



and substituting into the wave equation, Eq. ([3]) , we may 
group terms according to the order of the correction. To 



lowest order in the perturbation wc arrive at the stan- 
dard result that fields remain unchanged while the cigcn- 
frequencies shift according to 



Alo„i. = 





Ae\E,n) 


2 {E,„, 


e 


E„i) 



(8) 



For the applications that we have in mind the pertur- 
bation will typically be caused by chemical changes or 
reactions in the liquid, i.e. the perturbation will be non- 
zero in the liquid part of space only. Thus, in general we 
will encounter problems with Auj„i oc fm where 



{Em\£\Em) 



r, 



{Em\D,n) 



{EmlelEm) ^ (-Em -Dm)^ 



(9) 



is the filling fraction quantifying the relative optical over- 
lap with the liquid, i.e. < / < 1. In the second equal- 
ity of Eq. ^ we have introduced the displacement field 
|-D) = ej-E). The subscripts on the inner products em- 
phasize the spatial domain of integration in Eq. ([5]) . Note 
that compared to the definition normally used by the 
photonic crystal community [18] our definition quantifies 
the relative overlap with the liquid-part of space rather 
than the dielectric part given by 1 — / as easily seen from 
Eq. dH). 

For a fixed frequency the result of the perturbation 
may alternatively be viewed as a shift An in the Bloch 
wave vector k. Mathematically, the shift is calculated 
via the chain-rule 



AKr, 



dn 



Aujm = — (10) 

2 ei Vn,m 



diOm/ dn is the group velocity of the mth 



where Vg_ 
mode. 

In general, for enhanced light-matter interaction phe- 
nomena there is a balance between i) modifying the dis- 
persion (compared to propagation through the homoge- 
neous liquid) and ii) still maintaining a reasonable op- 
tical overlap with the liquid, i.e. keeping / of the order 
unity rather than close to zero. 



3 Dielectric properties of physiological liquids 

In this work we emphasize applications involving phys- 
iological liquids as can be found in biological systems. 
At low frequencies water is highly polarizablc because 
of the polar nature of the H2O molecules. In the static 
limit the relative dielectric function is around e 80. 
However, at optical frequencies the polar molecules move 
too slowly in respond to the temporally rapidly varying 
field and water becomes transparent to visible and in- 
frared radiation. The dielectric function is typically of 
the order e ^ (1.33)^ ~ 1.8. Compared to highly pu- 
rified distilled water, physiological liquids have natural 
contents of various salts in dissolved form. One example 
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could be potassium chloride (KCl) which has a high sol- 
ubility in water. The potassium (K+) and chlorine (Cl~) 
ions are highly mobile in water with a mobility around 
Mion ~ lO^^m^(Vs)^^. For ions carrying a charge Ze 
the ionic conductivity is given by 



a — ZeCion^J-ion 



(11) 



where Cion is the ionic concentration. Typically, the ionic 
concentration is of the order Qon ~ 1 mM yielding con- 
ductivities of the order Cion ^ 10~^S/m [45] . The elec- 
trical field will induce a current density J = aE in the 
liquid, but from an optical point of view the form of 
Eq. ([3]) is formally retained, but now with a complex 
valued dielectric function 



(12) 



Consequences of the imaginary part were discussed in 
Ref. [33]. We here emphasize the quite similar form of 
Eq. to the Drude model often employed for the re- 
sponse of metals at optical frequencies. For a liquid it 
will be an adequate description of its bulk properties, 
thus neglecting any surface chemistry as well as Debye- 
layer ion accumulation at the interfaces to other materi- 
als such as the high-index material forming the photonic 
crystal. Another issue could be electro-hydrodynamic 
phenomena where momentum is transferred from the 
ionic motion to the fluid (see e.g. Ref. [IS] and refer- 
ences therein), but at optical frequencies such phenom- 
ena are strongly suppressed, since the Debye screening 
layer forms too slowly in response to the rapidly vary- 
ing electrical field. It is common to introduce the Debye 
response time 

TD = lod^ = £il<y (13) 

and for typical electrolytes, tq is of the order a micro 
second corresponding to a Debye frequency in the mega- 
hertz regime. For optical frequencies, w ^ ujjj, it is thus 
fully adequate to treat the Ohmic damping due to the 
imaginary part in Eq. I|12p perturbatively. The damping 
is given by = 2Im{/«m} which from Eq. pO)) becomes 



a. 



(14) 



The attenuation will thus quite intuitively increase with 
a slowing down of the electromagnetic mode near photonic- 
band edges. Modes with a large optical overlap with the 
liquid of course suffer most from Ohmic dissipation as 
reflected by the proportionality of a to /. In a spatially 
homogenous liquid, Eq. (|14[) with / = 1 results in a 
damping of the order lod/c corresponding to a typical 
decay length of the order 1 m. For applications of liquid- 
infiltrated photonic crystals in lab-on-a-chip applications 
it thus seems adequate to treat typical physiological liq- 
uids and the high-index dielectric material at an equal 
footing, i.e. both materials are to a first approximation 
assumed lossless. 




Fig. 3 TM modes for a square-lattice of pillars with di- 
ameter d/A = 0.5 and dielectric function £d — 10.5. Panel 
(a) shows results for the non-infiltrated photonic crystal and 
panel (b) is for the same photonic crystal infiltrated by a 
hquid with dielectric function e; = (1.33)^. 



4 Classification of photonic bands 

The majority of experimental work on liquid-infiltrated 
photonic crystals has so far emphasized void-like struc- 
tures. The pillar-like structures may be viewed as the ge- 
ometrical inverse with the high and low-index materials 
being swapped. From Babinet's principle [47| one might 




Fig. 4 TE modes for a triangular-lattice of voids with diam- 
eter d/A = 0.96 in a dielectric material with Ed = 10.5. Panel 
(a) shows results for the non-infiltrated photonic crystal and 
panel (b) is for the same photonic crystal infiltration by a 
liquid with dielectric function ei = (1.33)^. 
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imagine the two types of photonic crystals to have some- 
what related properties. However, the source of forma- 
tion of photonic band gaps are quite different in nature; 
for void-like structures the band-gaps arise from multiple 
reflections, i.e. Bragg scattering while band-gaps relate 
to Mie resonances for pillar- like structures |48ll49j . The 
sensitivity to infiltration by liquids is also somewhat dif- 
ferent in terms of polarization. The discussion of dielec- 
tric versus air bands in Ref. [18| is naturally extended to 
the present context where we will refer to the bands as 
dielectric or liquid bands. Fig. [3] illustrates the typical 
TM properties of a square-lattice pillar-like structure. 
The band structure has been obtained with the aid of a 
plane- wave simulation [44j . The lowest band is obviously 
a dielectric band with / ~ as clearly seen by the mod- 
est change upon infiltration. The second band is on the 
contrary a liquid band and the higher value of / causes a 
pronounced frequency shift of the upper band-gap edge 
upon infiltration by the liquid. Due to the reduced in- 
dex contrast, especially high-order band-gaps may close 
upon liquid infiltration as seen by comparing panels (a) 
and (b) in Fig. E] 

Typical TE properties of a triangular-lattice void- 
like structure is illustrated in Fig. |4l The band structure 
has been obtained with the aid of a plane-wave simu- 
lation [33]. Again, the lowest band is predominantly a 
dielectric band with / ~ 0.4 as clearly seen by the mod- 
est change upon infiltration, though this is not as pro- 
nounced as for the pillar-like structure in Fig. O Like- 
wise, the second band is a liquid band with a high value 
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Fig. 6 Slow-light enhanced absorbance. Panel (a) shows the 
band structure for propagation of TM polarized light along 
the rX direction in a line-defect waveguide in a square lat- 
tice of period A with dielectric rods of diameter d — OAA 
and e = 10.5 with the structure infiltrated by a liquid with 
£ = (1.33)^. The complete photonic band gap of the photonic 
crystal is indicated by yellow shading while grey shading in- 
dicates the finite density-of-states in the photonic crystal due 
to the projected bands in the Brillouin zone. Panel (b) shows 
the corresponding enhancement factor 7 which exceeds unity 
over the entire bandwidth. The right y-axis shows the results 
in terms of the free-space wavelength when results are scaled 
to a structure with A = 420 nm. 
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Fig. 5 Filling fraction for first dielectric (red line) and liquid 
(blue line) bands. Panel (a) shows results corresponding to 
the band structure in panel (b) of Fig.|3]while panel (b) shows 
results corresponding to panel (b) of Fig. |4] The insets show 
the electrical field intensity \E\'^ at various high-symmetry 
points for the two type of bands. 



of / causing a pronounced shift of the upper band-gap 
edge upon infiltration by the liquid. 

Modes with a liquid-band behavior are of course of 
primary interest since the high electrical field intensity 
in the liquid allows for optofluidic tuning of the disper- 
sion or likewise highly sensitive optical monitoring of 
chemical changes in the fluid. Fig. [5] shows results for 
the filling fraction / along with field-intensity plots for 
the dielectric and liquid bands associated with the band 
structures in panels (b) of Figs. [3] and [H High values of 
/ correlates with a pronounced field localization in the 
liquid while low values is associated with a marked field 
overlap with the dielectric. 



5 Applications in sensing and quantitative 
analysis 

5.1 Beer- Lambert- Bouguer absorption 

The Beer-Lambert-Bouguer absorbance technique was 
developed two centuries ago and has since then become 
one of the classical optical workhorses in analytical chem- 
istry [11] . The basic principle of an absorbance measure- 
ment is illustrated in Fig. 1. The optical probe, with 
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Fig. 7 Transmission of TM polarized light incident on a liquid-infiltrated photonic crystal based on a square-lattice of pillars 
with diameter d/A — 0.4 and dielectric function Sd ~ 10.5. Panel (a) shows the spectral transmission for the non-infiltrated 
photonic crystal (dashed line) and when infiltrated with a liquid with n; = 1.33 (solid green line) and n; = 1.5 (solid red line), 
respectively. Panel (b) shows the spectral shift of the upper (solid line) and lower (dashed line) transmission band edges for 
varying values of the liquid refractive index n; . 



intensity Iq, is incident on a liquid sample with absorp- 
tion parameter ai reflecting a certain concentration C 
of some chemical species to be quantified. Typically the 
chemicals will be dissolved in a liquid, but for certain 
applications gas and even transparent solid phases are 
in principle also possible. To lowest order in the concen- 
tration we have 

ai^^ViCi, (15) 

i 

with the Ci being the concentration of the absorbing 
bio-molcculcs of the ith type with molar extinction co- 
efficient rji. 

If we neglect coupling issues, the transmitted inten- 
sity / will then, quite intuitively, be exponentially damped, 

I = lo exp(-7a;L). (16) 

We have, by hand, introduced 7 as a dimensionless mea- 
sure of possible dispersion enhanced light-matter inter- 
actions. For a uniform medium we of course have 7 = 1 
and the expression then often goes under the name of 
Beer's law. For a dilute solution ai correlates linearly 
with the concentration of the absorbing chemical species 
and Beer's law then provides simple optical means for 
detecting and quantifying the concentration of chemi- 
cal solutions [11]. Evidently, the effect relies heavily on 
having a sufficiently long optical path length and the 
longer L is the lower a concentration can be monitored 
for a given sensitivity of the optical equipment measur- 
ing I/Iq. Quite obviously, lab-on-a-chip implementations 
of chemical absorbance cells are facing major challenges 
due to the miniaturization as exemplified by recent ex- 
perimental work by Mogensen et al. [15j . Here, we out- 
line our recent suggestion to use slow-light phenomena 



to compensate for the reduced optical path length [35] , 
In order to rigorously derive the enhancement factor 7 
for a liquid-infiltrated photonic crystal we apply per- 
turbation theory, Eq. (|10|1 . to the problem of a small 
imaginary part iAsi in the liquid. Comparing the cor- 
responding damping cxm = 2Im{K„i} to the damping in 
a spatially homogenous liquid, obtained from Eq. (jlOp 
with / = 1 and Vg = c/ni, we arrive at an enhancement 
factor 7„; = a^n/ai given by 

7m = /™x^. (17) 

The interpretation of this result is straightforward; slow- 
light propagation, Vg <C c, will significantly enhance the 
absorbance. The presence of the filling factor / is also 
quite obvious since only the relative fraction of the in- 
tensity residing in the liquid will contribute to the ab- 
sorption. Finally, we note that for a spatially homoge- 
neous liquid we get 7 = 1 as required by the definition 
in Eq. p6p . In general, the enhancement factor has a 
lower bound, 7 > 0, while there is in principal no upper 
bound since the inverse group velocity diverges for Bloch 
wave vectors near the Brillouin zone (such as the X or M 
points) or near the center of the zone (the F point), see 
Fig- El In Fig. [5] we show results for a line-defect waveg- 
uide in a square lattice of period A with dielectric rods 
of diameter d = OAA and e = 10.5. As expected, the 
pillar-like structure supports a mode with a very high 
optical overlap with the liquid, / ~ 1, combined with a 
reduced group velocity over the entire band width of the 
waveguide transmission window, see panel (a). Further- 
more, at the r point the inverse group-velocity diverges 
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which is also the case near the upper waveguide trans- 
mission edge. In fact, the enhancement factor is of the 
order or larger than 2 over the entire bandwidth. En- 
hancement by an order of magnitude is easily achieved 
near the transmission band edges. 



5.2 Refractometry 

Refractometry in its basic form is another widely used 
method which most often is a qualitative measurement 
indicating chemical changes contrary to e.g. absorbance 
measurements which are of quantitative nature. Chem- 
ical reactions and changes most often cause a modifica- 
tion Am of the refractive index, thus altering the op- 
tical response of e.g. a reflection or transmission mea- 
surement. Chemical changes may thus be monitored by 
watching the spectral shifts of e.g. transmission reso- 
nances j33l[36] . Obviously, the larger a spectral shift the 
more sensitive is the method. As for the previous prob- 
lem we will apply perturbation theory to study the spec- 
tral sensitivity due to a small perturbation Aei ~ 2niAni. 
From Eq. ^ we get 



AuJr, 



— .fm X 



Ae^ 

2e, 



(18) 



or introducing the free-space wavelength A = l-Kcjio we 
may express the sensitivity in the form of a wavelength 
shift A\ due to a shift in liquid refractive index Ani^ 



AXm 

Am 



fm X 



A. 



ni 



(19) 



The interpretation is again quite intuitive. In a homoge- 
neous liquid we have a dispersion 1/A oc uj{k) = {c/m)k 
so that AX/ Am = X/m which is obviously the ulti- 
mate sensitivity one can have. However, transmission 
through homogenous media is not associated with any 
pronounced resonances or cut-offs. Transmission exper- 
iments reflect the photonic dcnsity-of-states and strong 
resonances may be achieved by infiltrating a photonic 
crystal with the liquid [34] . In that case the sensitivity is 
naturally corrected by the filling factor / as shown by the 
rigorous perturbative analysis in Eq. (fT9|) . From Eq. (fT9|) 
the working wavelength is of course important to the 
sensitivity. However, this is a somewhat trivial effect for 
refractometry in general, while the dependence on fill- 
ing fraction in Eq. (fT9|) is really the limiting factor. The 
performance of different structures, ranging from evanes- 
cent field sensing devices to liquid-infiltrated photonic 
crystals, may thus be compared by comparing the cor- 
responding filling fractions. Such relative comparisons 
has previously been suggested in the literature [30| . but 
Eq. (fT9|) does in fact allow for an absolute definition of 
the sensitivity. 

Evanescent field-sensing devices are characterized by 
an almost vanishing / which will typically only be of the 
order a few percent. Depending on the design, void-like 



photonic crystals offer moderate-to-high filling fractions 
depending on the void-to-solid volume ratio, and simi- 
larly for pillar-like photonic crystals. 

Figure [7] shows a proposal of refractometry with a 
liquid-infiltrated photonic crystal. The simple structure 
consists of a square-lattice of pillars with diameter d/ A = 
0.4 and dielectric function Sd = 10.5. Transmission spec- 
tra for TM polarized light have been obtained by means 
of a two-dimensional finite-difference time-domain (FDTD) 
method, see Ref. [33] and references therein. Panel (a) 
shows the spectral transmission for the non-infiltrated 
photonic crystal (dashed line) and when infiltrated with 
a liquid with m ~ 1-33 (solid green line) and m — 1-5 
(solid red line) , respectively. From the band structure in 
Fig. [3] one would expect a suppression of the transmis- 
sion for frequencies inside the photonic band gap which 
correlates nicely with the observed transmission. Fur- 
thermore, when infiltrating the photonic crystal with liq- 
uid, the spectral position of the upper transmission band 
edge is expected to be more sensitive to the perturbation 
than the lower transmission band edge. Panel (b) shows 
the spectral shift of the upper and lower transmission 
band edges for varying values of the liquid refractive in- 
dex m ■ The upper curve does indeed have a larger slope 
than the lower curve which according to Eq. (|19p sup- 
ports the picture of transmission through a liquid band 
and dielectric band, respectively. The transmission re- 
sults agrees quantitatively with the band-structure the- 
ory as studied in detail in Ref. [33] and from Eq. (fT9|) 
we estimate that / ~ 80% for the transmission through 
the liquid band. 

In Fig. [8] we show another implementation of re- 
fractometry employing a liquid-infiltrated photonic crys- 
tals wave-guide structure based on a triangular-lattice 
of voids with diameter d/A = 0.6 and dielectric func- 
tion Ed = 10.5. Selective filling of the defect wave-guide 
structure is of course crucial, but nanofluidic tuning has 
recently been demonstrated experimentally in a similar 
structure [3|. Panel (a) shows the spectral transmission 
of TE polarized light for the non-infiltrated photonic 
crystal (dashed line) and when infiltrated with a liquid 
with m — 1-33 (solid green line) and n/ — 1.5 (solid red 
line), respectively. Panel (b) shows the spectral shift of 
the upper and lower transmission band edges for varying 
values of the liquid refractive index n;. Again, the upper 
curve corresponds to transmission through a band of pre- 
dominantly liquid-band nature with / ~ 40% estimated 
from Eq. (fTO]) . 

Above we have emphasized transmission properties, 
but spectral shifts could alternatively be monitored in a 
reflection configuration. 

5.3 Cavity based sensing 

Obviously, cavity modes may be used for refractometry 
where Eq. also applies. Defect structures in pho- 
tonic crystals are quite interesting in this context and 
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Fig. 8 Transmission of TE polarized light incident on a liquid-infiltrated photonic crystals wave-guide structure based on a 
triangular-lattice of voids with diameter d/ A — 0.6 and dielectric function Ed ~ 10.5. Panel (a) shows the spectral transmission 
for the non-infiltrated photonic crystal (dashed line) and when infiltrated with a liquid with ni — 1.33 (solid green line) and 
ni = 1.5 (solid red line), respectively. Panel (b) shows the spectral shift of the upper (solid line) and lower (dashed line) 
transmission band edges for varying values of the liquid refractive index n; . 



in recent experiments Chow et al. have pursued this di- 
rection |23| . For refractometric applications the quality 
factor Q of a resonance should be sufficiently high to 
easily monitor and resolve even small spectral shifts of 
the resonance. Here, we would like to emphasize another 
potential usage of high-Q resonators somewhat resem- 
bling the Beer-Lambert-Bouguer absorption technique. 
For the slow-light enhancement we note that the diver- 
gence of the inverse group velocity near the edge of the 
Brillouin zone is basically a standing wave phenomenon. 
Thus, cavity modes may in a similar way be used to en- 
hance the light-matter interaction time r, but contrary 
to the case of Eq. ^ the Wigner-Smith delay time now 
relates to the cavity quality factor Q. In the following 
we consider an isolated Lorentzian resonance. 



p,n{i^) = --Im 

TT 



1 



(20) 



with the line width corresponding to an intrinsic 
quality factor ~ LUm/rm- We imagine that the in- 
trinsic quality factor includes everything but the broad- 
ening due to absorption caused by the presence of bio- 
molecules in the liquid. The total quality factor Q may 
then be written as 

I Qahs 7 Qahs Qin^ 

Q = {Qrn+Q:^Ly ^\ (21) 

[ Qin , Oil ^ 0, 

where Qabs is the contribution from absorption. High-Q 
cavities would be of prime interest to sensing since the 
line width would be dominated by the presence of ab- 
sorbing bio-molecules rather than the intrinsic line width 
associated with the localization of light. Obviously, the 



higher a Qin the lower a concentration of molecules can 
be detected. In order to quantify this further we con- 
sider the limit of weak absorption in the liquid and ap- 
ply perturbation theory. Calculating the imaginary shift 
of uJm due to Ae = iAei from Eq. ([8]) and substitut- 
ing into Eq. (j20p we identify an additional broadening 
corresponding to 

27rni 



/abs 



(22) 



where we have introduced the absorption parameter a; 
of the liquid. The form of Eq. resembles the form de- 
rived for e.g. homogeneous micro-sphere resonators |50| 
except for the presence of the filling factor /. Obvi- 
ously, the lower an overlap with the absorbing liquid 
the smaller an additional broadening of the intrinsic res- 
onance line width. 

If we consider just one type of molecules, see Eq. 
the criterion in Eq. pT|) means that the lowest detectable 
concentration is bound by 

1 



C > f 



(23) 



with 77 being the molar extinction coefficient. As an ex- 
ample, for a typical value rj lM~^/^m~^ a detection 
in the micro- molar regime would require / x Qm ^ 10^. 

Cavity designs optimizing the product / x are ob- 
viously needed, but it is also well-known that for slab-like 
photonic crystals the typical high-Q cavities will have 
a quite modest / to prevent degradation of Q by out- 
of-planc leakage [5T|[52] though some designs may get 
around this [551l54j . Indeed, the refractometric measure- 
ments in Ref. suggest that / ~ 20% for their cav- 
ity mode with Q ~ 400. Cavities in three-dimensional 
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photonic crystals with photonic-band gap confinement 
in all three directions could be an interesting scenario. 
As an alternative, recent experiments with liquid micro- 
droplets facilitate whispering-gallery modes with 
/ ~ 1 and potentially the intrinsic Q-factor can be made 
very high due to the absence of surface roughness at the 
almost perfect liquid-air interface which is smooth on 
the molecular length scale. Ref. [56] reported Q ^ 10^ 
in a fiber-coupled experiment. Experimentally one might 
use the above scheme to quantify concentrations e.g. by 
time resolved ring-down spectroscopy or by other means 
such as transmission experiments typically employed for 
characterization of high-Q resonators [57] . 

6 Conclusion 

With the increasing emphasis on miniaturization of chem- 
ical analysis systems there is presently a strong effort 
devoted to integrating microfluidics and optics in lab- 
on-a-chip microsystems. We have outlined possible ap- 
plications of optofluidics in this context and used sim- 
ple theoretical approaches to illustrate how light-matter 
interactions may be enhanced in liquid-infiltrated pho- 
tonic crystals supporting highly dispersive modes and 
pronounced resonances. 

Our formalism combines electromagnetic perturba- 
tion theory with full-wave electromagnetic simulations. 
The latter is used in calculating the dispersion proper- 
ties of ideal (lossless) liquid-infiltrated photonic crystals 
while the former is used systematically to predict the 
spectral effects and changes caused by chemical pertur- 
bations of the infiltrated liquid. 

In this paper we implicitly focused on near-infrared 
wavelengths with values for the high-index material re- 
sembling silicon. However, we emphasize that the general 
concepts apply equally well to other wavelength regimes 
such as the visible regime as in Ref. [10] or even the ter- 
ahertz regime emphasized in Refs. [26 l [30 l l3T| . The phe- 
nomena that we have discussed do to some extend rely 
on a sufficient index contrast between the liquid and the 
high-index materials. While polymers are certainly at- 
tractive materials for future low-cost disposable lab-on- 
a-chip devices they might not provide a sufficient index 
contrast for all of the applications that we have discussed 
and the choice of wavelength and materials should of 
course reflect the particular applications in mind. 

With this theoretical paper we have sought to empha- 
size the potential benefits from using liquid-infiltrated 
photonic crystals for sensing purposes in lab-on-a-chip 
systems such as label- free sensing of bio- molecules. While 
photonic crystal concepts are generally considered a pro- 
mising route to highly sensitive devices, our theory clearly 
demonstrates an inevitable challenge; the photonic crys- 
tal design needs to balance between having a strong opti- 
cal proximity from the photonic crystal structure and at 
the same time supporting a decent optical overlap with 



the liquid. On the other hand, the huge design space of- 
fered by photonic crystals seems to allow us to explore 
and cope with this challenge. In conclusion we believe 
that photonic crystals may be used to facilitate enhanced 
light-matter interactions which may be traded for yet 
smaller miniaturized systems or for increased sensitivity 
of existing devices. 
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